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Abstract 

Log-likelihood ratio (LLR) computation for non-binary modulations over fading channels is compli- 
cated. A measure of LLR accuracy on asymmetric binary channels is introduced to facilitate good LLR 
approximations for non-binary modulations. Considering piecewise linear LLR approximations, we prove 
convexity of optimizing the coefficients according to this measure. For the optimized approximate LLRs, 
we report negligible performance losses compared to true LLRs. 
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I. Introduction 

It is well known that soft-decision decoding algorithms outperform hard-decision decoding 
algorithms. In soft-decision decoding, reliability metrics are calculated at the receiver based on the 
channel output. The decoder uses these reliability measures to gain knowledge of the transmitted 
codewords. The superiority of soft decoding comes at the expense of higher complexity. 

Log-likelihood ratios (LLRs) have been shown to be very efficient metrics for soft decoding of 
many powerful codes such as the convolutional codes lU, turbo codes [2J, low-density parity-check 
(LDPC) codes [|3l. LLRs offer practical advantages such as numerical stability and simplification 
of many decoding algorithms. Moreover, due to some properties of the probability density function 
(pdf) of the LLRs, such as symmetry and invariance, LLRs are used as convenient tools for the 
performance analysis of binary linear codes Il4l-|l6l. On many communication channels, even 
for binary modulations, i.e., binary phase-shift keying (BPSK), channel LLRs are complicated 
functions of the channel output |I71. This fact greatly increases the complexity of the LLR 
calculation modules in the decoder causing decoding delays and power dissipation. In high speed 
wireless transmissions, the decoder may not be able to handle this complexity. Thus, for an efficient 
implementation of the decoder, approximate LLRs should be considered. 

Approximate LLRs have been previously used in the literature [TTII- lfTOll . Piecewise linear LLRs 
have been suggested in [|9l for soft Viterbi decoding of convolutional codes in the HIPERLAN/2 
standard [UJ. The presented method uses the log-sum approximation which is quite accurate 
at high signal-to-noise ratio (SNR). Moreover, it assumes that perfect channel state information 
(CSI) is available at the receiver. In [|12|. linear LLRs have been used for BPSK modulation on 
uncorrelated fading channels without CSI and a measure of LLR accuracy has been introduced. 
Using that measure, linear LLR approximating functions have been designed with almost no per- 
formance gap to that of true LLR calculation. The proposed measure, however, is only applicable 
to symmetric channels and BPSK. With non-binary modulations, used in most practical systems, 
a linear approximation of bit LLRs is not always possible. Moreover, the equivalent bit-channels 
are asymmetric. That is, the LLR accuracy measure of [12] is not applicable. 

In this work, we seek approximate LLRs for non-binary signalling over uncorrelated fading 
channels. Thus, we need to generalize the accuracy measure of [12] to asymmetric channels. 



Moreover, we must consider non-linear LLR approximating functions. To compute LLRs at bit 
levels individually, we consider bit-interleaved coded modulation (BICM) lfT3l . BICM is a well- 
known bandwidth efficient scheme for binary codes on fading channels. 

While our approaches are general, to demonstrate our methods, we focus on piecewise linear 
approximations. Such approximations are easy to implement and we observe that, when the 
parameters are optimized, they perform close to true LLRs. We prove that the optimization of 
these piecewise linear approximations is convex, thus, very efficient. Due to the close-to-capacity 
performance of LDPC codes on many channels (HI, UM, lESIl, we employ LDPC-coded BICM |[T6l 
to show that even with the proposed approximation, close-to-capacity performance is obtained. 

IL Problem Definition 
Consider a flat slow-fading environment where the received signal is expressed as 

y = r ■ X + z, 

where x is the complex transmitted signal chosen from the signal set A:' C C of size l^:*! = 2"^, 
r > is the channel fading gain with arbitrary pdf p{r), and z is the additive noise which is a 
complex zero-mean white Gaussian random variable with variance 2cr^. 

A. LLRs for equivalent bit-channels 

Using the BICM scheme [fT3l . the information sequence is first encoded by a binary code. Next, 
the coded sequence is bit interleaved and is broken into m-bit sequences which are then Gray 
labeled onto signals in X and transmitted on the channel. Assuming ideal interleaving, the system 
can be seen equivalently as m parallel independent and memoryless binary-input bit-channels. In 
the receiver, based on y, LLRs are computed for each bit-channel independently from other bits. 
These LLRs are then de-interleaved and passed to the decoder. 

When the channel fading gain r is known at the receiver for each channel use, the true LLR 
for the ith bit-channel, assuming uniform input distribution, is calculated as 

P(y|6'(x) = l,r) ExG;t^Hy|x,r) 
where i E {1, • • • ,rn}, 6*(x) is the ith bit of the label of x, is the subset of signals in X 
where b*(x) = w, and the conditional distributions are given by p(y|x, r) = j^exp (— ) 



or exp ( — -^^2^1^) when the signal set is real. Also, gr\y) represents L^^^ as a function of 
y when r is known. When r is not known at the receiver, the true LLR is calculated as 

i... . ^W^-W - °) . log P^J^^^ ^ ,<..(y). (2, 

P(y|b'(x) = l) Exe;t^Hy|x) 

where p(y|x) = Jp°° 2;^^^? (— ^^2^2^^ )p(r)dr, and fi'*-*-'(y) represents L*^*) as a function of y. 

As can be seen from ([T]) and both (7r*''(y) and (y) are usually complicated functions 
of y. Thus, approximate LLRs (^^*^(y) and g^^^y)) are of practical interest. One approximation 
which is useful at high SNR is obtained by the log-sum approximation: log J2k ~ maxjt log zt- 
This approximation is good when the sum is dominated by a single large term. Thus, 

max„p v» p(y|x, r) 

9Hy) = ^og — ^"'^°7 , (3) 

^ 1 ^^^yceXiPiyl^) , max^g v. p(y |x, r)p(r)dr 

y = log ° I . = log \oo . I TTTT- 

max^g;t'^P(y|x) max^g;,., Jq p(y|x, r)p(r)dr 

The log-sum approximation is particularly useful when CSI is available at the receiver, where 

^ leads to piecewise linear LLRs which can be efficiently implemented [l9l|. However, with no 

CSI, the log-sum approximation is no longer piecewise linear and in fact involves complicated 

integrations. The focus of our work is on cases that CSI is not available. Nonetheless, we seek 

approximate LLRs which are piecewise linear functions of y. 

Now, consider a general LLR approximation function parameterized by set of parameters Ai 

This function maps the complex received signal y to real-valued approximate LLR for the ith bit- 
channel. Clearly, it is desired to choose Ai such that accurate LLR estimates are obtained. To this 
end, we first extend the LLR accuracy measure introduced in lfT2l to non-symmetric bit-channels 
obtained via BICM. We then use this measure for optimizing the LLR approximating functions. 

B. LLR accuracy measure for asymmetric binary-input channels 

In this section, we generalize the LLR accuracy measure of lfT2]| to asymmetric channels. To this 
end, we consider the pdfs of the ith bit-channel LLR conditioned on the transmitted bit h E {0, 1}, 
defined as p^f^\l) = E^(^p^i[p{L^^^ = /|x)]. Using these LLR pdfs, it is possible to calculate the 
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capacity of each asymmetric bit-channel and thus the BICM scheme. By capacity, we mean the 
mutual information between the input and output of each bit-channel when its input b is equally 
likely or 1. The capacity of the ith bit-channel is given by [fTTll 

a=l-^l log2(l + e-')rf^(/)d/ - i 1 log2(l + e^)p^;\l)dl. (4) 

Thus, the capacity of the BICM is found by C = Xll^i 

When instead of the true LLRs, approximate LLRs L^^^ are used, their pdfs are given by 
pI^\1) = E-^^;:^;i[p{L^^^ = /|x)] for b G {0, 1}. Inserting these approximate LLR pdfs in (HI) we get 

^ = E^^ = E(i-^/ + e-^)p^^\l)dl -\j \og,{l + e')pfmi) ■ (5) 

The following theorem, proved in Appendix A, shows C can be used as an LLR accuracy measure. 

Theorem 1. The maximum of C in ((5]) is equal to C which is achieved by true LLRs (no LLR 
approximation can result in C > C). 

Notice that for a symmetric channel (i.e., Pi\l) = e~'po^(/)), ([5]) reduces to the LLR accuracy 
measure of [fT2|. Using similar arguments, it can be shown that as the approximate LLRs drift 
away more from the true LLRs, AC = C — C gets larger. Thus, (7 is a measure of the accuracy of 
the approximate bit LLRs. Since bit-channels are independent, we maximize each Ci individually. 
Thus, for each bit-channel i, assuming a class of approximating functions fi'^^(y), we find: 

.4°''* = argmaxQ, (6) 
s.t. *,(A)=o 

where ^i{Ai) = denotes the constraints imposed on Ai (e.g., to preserve continuity). 

C. LLR approximating functions 

It is evident from ([T]) and © that true LLR functions depend on the signal set X. Thus, choosing 
appropriate class of approximating functions also depends on A*. In this letter, we consider non- 
binary AM and rectangular QAM. Our optimization approach, however, is general. Moreover, we 
put our focus on piecewise linear approximate LLRs. Clearly, the optimization problem Q can be 
solved for any other approximating function. Using piecewise linear approximations has benefits 
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such as simplicity of demodulator and (as will be shown) convexity of the optimization problem. 
Numerical results verify that the obtained performance is also very close to true LLRs. 

By viewing a complex variable as a two-dimensional vector y = (Re{y}, Im{y}), a piecewise 
linear function of a complex variable is defined as follows. First, the complex domain C is divided 
into a finite number of regions Ci, C2, . . . , Cat by a finite number of one-dimensional boundaries. 
Then the function is represented by /(y) = (a^, y) + bk for any y G Ck, where (, ) denotes the 
inner product of two-dimensional vectors, i.e., (afc,y) = Re{afc}Re{y} + Im{afc}Im{y}. Thus, 

k=l 

where A^*^*^ is the number of segments of the piecewise linear function. At is the set of all a^*\ 
b^i!:\ and 1(.) is the indicator function. The parameters are chosen to preserve continuity over y. 

The following theorem, proved in Appendix B, indicates that optimizing a piecewise linear 
approximating function according to Q is a convex optimization problem. 

Theorem 2. Assuming that approximate bit LLRs are calculated by ([7|) for i = 1, . . . ,m, and 
assuming fixed C^*'* for k = 1, . . . , N^^\ Ci is a concave function of a^^ and h^^ for all k. 

Using (|7]), ^ can be numerically solved as follows. For a given SNR, and assuming fixed 
C^*^'s, Ci can be computed by first computing Pq\1) and Pi\l) for given a^*'''s and h^l^'s, and 
inserting them in dS]). Since Ci is a concave function of a^*^'s and l>{^'^ ^'^^ '^he constraints 
are linear, maximizing Ci can be done efficiently using numerical optimization techniques. The 
proper number of regions A^*^*^ is based on the affordable complexity and the curve of true LLRs. 
Optimizing the regions can be done through search. As will be seen in the next section, usually 
the size of the parameter sets is small and due to the symmetry in the LLR calculation, many 
parameters are equal to each other which can further reduce the number of unknown parameters. 

in. Proposed Approach and Examples 

Now, we describe the proposed method through examples of real and complex signal constel- 
lations. 
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Example 1: Consider 8-AM constellation with Gray labeling shown in Fig.[T]on the normalized 
Rician fading channel. On this channel, p{r) = 2r{K+l)e-^^^^^-^'^>^ha{2r^K{K + 1)), where 
K is the Rician K-factor, and /o(-) is the zero-order modified Bessel function of the first kind. 

Using ([U) and true LLRs are calculated for i = 1,2, 3, and they have been plotted versus 
y for extreme values of K in Fig. [2l Considering the general model of (|7]), and the symmetry in 
the true LLR functions, we propose the following piecewise linear LLR approximations 

L^''^=9^Il{y) = {a?y + b?%y<o) + {a^2^y + b?%o<y) = -a?\y\ + b?, (9) 

L^'' = a'lliy) = ia?y + bf%y,,^^, + {afy + b^')l^^,^,^ 

+ (af y + &f )l(o<,<43)) + (af 2/ + )l(ef)<,), (10) 

where due to the symmetry of the LLRs, we have assumed in ^ that a^2^ = —af^ and b^^'^ = bf\ 
Also, in (Uni), we have af = -af\ = -af , bf^ = bf , bf = bf, and = = 
-cf = -cf . Thus, Ai = {aS'^}, A2 = {af\bf^}, and ^3 = {af\ a^^\bf ,bf\ cf^}. It is 
evident that symmetry reduces the number of unknown parameters. We also impose $(^3) = 
c[ \af^ — ttg'^'') + bf^ — b^2 ^ = to preserve continuity in (flOl) . Fig. [21 shows that these piecewise 
linear functions better approximate the true LLRs when K increases. 

For a given SNR, we optimize the parameter sets Ai, A2, and ^3, by solving Numerical 
results confirm that Cuiax = J^iLi max^i Ci is always very close to the capacity of BICM 
employing true LLRs, i.e., C. For example, for K = 0, Cmax = 0.851 and C = 0.855 bits 
per channel use at SNR= 5.00 dB, and Cmax = 1-544 and C = 1.553 bits per channel use at 
SNR= 30.00 dB. When K increases, AC = C — Cmax becomes even smaller. 

To evaluate the decoding performance of the optimized piecewise linear approximations, we 
compare the decoding threshold of LDPC codes and their bit error rate (BER) under approximate 
and true LLRs. The decoding threshold can be found by density evolution 051, lfT4l and by using 
the technique of i.i.d. channel adapters [fT6l which provides the required symmetry conditions. 

As an example, consider (3, 4) -regular LDPC codes on the normalized Rayleigh channel (equiv- 
alent to K = G) with 8-AM signalling of Fig. [T] By using the approximating functions of (l8])-(fT0l). 
we find the decoding threshold of the code under optimized LLR parameters reported in Table HI 
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The decoding threshold given by density evolution is 7.88 dB while under true LLR calculation 
of dll), the decoding threshold is 7.85 dB showing only a 0.03 dB performance gap. 

To see how the piecewise linear LLR calculation affects the BER performance, we simulate a 
given LDPC code on the normalized Rayleigh fading channel. In Fig. [3l the performance of a 
randomly constructed (3, 4)-regular LDPC code of length 15000 is depicted in two cases: once 
decoded using true LLRs of and once with the piecewise linear approximation of (f8l)-(fT0l) and 
the optimized parameters reported in Table HI It should be noted that the parameters are optimized 
once at the decoding threshold and are kept fixed at the receiver for other SNRs. It is seen that the 
performance of the optimized approximate LLRs is almost identical to that of the more complex 
true LLRs although parameters are only optimized for the worst SNR. 

Also, to show that optimizing C is meaningful in terms of the maximum transmission rate 
achievable by the piecewise linear LLRs, we optimize the degree distributions ifTSl of LDPC 
codes under our approximate LLRs. At SNR = 21.02 dB, the capacity of BICM under true LLRs 
is C = 1.500 in the absence of CSI when K = 0. Since m = 3, then the maximum binary 
code rate achievable on this channel is 0.500. At this SNR, solving ^ gives Cmax = 1-493 and 
the parameters reported in Table H Now, by using the designed piecewise linear approximation, 
assuming a fixed check node degree of 8 and maximum variable node degree of 30, an irregular 
LDPC code is designed. The variable node degree distribution of this code is X{x) = 0.250x + 
0.217x2+0.221x^+0.048x^+0.119x22+0.145x29, and the code rate is R = 0.490. Thus, the proposed 
approximate LLRs can achieve rates very close to the capacity of BICM under true LLRs. 

Example 2: Now consider a 16-QAM constellation with Gray labeling as depicted in Fig. [T] 
Using the general piecewise linear model of (|7]), due to the symmetry and the similarity of the 
bit LLR functions, we propose the following LLR approximations: 

i'"= 9A,{y) ='4"Re{y}, (11) 

k=l 

= Re{af)}|Re{y}| + Im{af ^}|Im{y}| + bf\ (12) 

(2) (2) 

where , . . . , C4 are the four quadrants of the complex plane. It should be noted that bit- 
channel LLR calculations are similar for bit 1 and 3, and for bit 2 and 4 except that the real 
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and imaginary parts of y are swapped, i.e., L*^^) = a^^^Imjy} and L^^^ = Re{a^^''}|Im{y}| + 

(2,) (2f) fl) (2f) (2f) 

lm{a\ "^jlRejy}! +b\' . Thus, it is enough to optimize Ai = {a\ '} and A2 = {8l\ ,b\ '}. 

Again numerical results suggest that the gap between Crnax and the true BICM capacity is 
always small. For example, when = 0, we have Cmax = 1-074 and C = 1.097 bits per channel 
use at SNR=3.00 dB. Again, the gap becomes smaller when K increases. 

To investigate the performance of 16-QAM signalling under the proposed approximate LLRs, 
we compare the decoding threshold and BER of (3, 4)-regular LDPC codes on the Rayleigh 
fading channel under true and approximate LLRs. Density evolution gives a decoding threshold 
of 5.02 dB under approximate LLRs with the optimized parameters of Table HI Under true LLR 
calculation, the decoding threshold is 4.83 dB. As a result, approximate LLRs show about 0.19 dB 
performance gap to true LLRs. The BER comparison is depicted in Fig. [3l It is worth mentioning 
that this gap can be further reduced by proposing piecewise linear LLRs with more segments. 

IV. Conclusion 

LLR computation for equivalent bit-channels of a non-binary modulation is generally compli- 
cated. On fading channels, when the channel gain is unknown, this problem is further intensified. 
Noticing that the equivalent bit channels were asymmetric, in order to find good approximate LLRs, 
we proposed an LLR accuracy measure for binary asymmetric channels. This accuracy measure 
can be used to optimize the parameters of any approximating function. We used our accuracy 
measure to optimize piecewise linear LLR approximations. By using LDPC-coded BICM, we 
showed that the performance loss under the optimized piecewise linear approximation was very 
small. We also showed that under approximate LLRs, asymptotic irregular LDPC codes having 
rates very close to the capacity of BICM under true LLRs can be obtained. Our solution can also 
be applied to other coding schemes which use LLRs such as the convolutional and turbo codes. 

Appendix A 
Proof of Theorem [H 

Consider an arbitrary discrete binary-input memoryless channel whose output alphabet is non- 
binary. The channel input is 2; G {0, 1}, and its output is y G {yj\l < i < M}. Let us define 
P{yj\x = 0) = pj and P{yj\x = 1) = qj where J2jLiPj = J2jLi1j = 1- The true LLR value. 
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when y = yj is observed at the channel output and the binary inputs are equiprobable, is 



(13) 



Thus, the true LLR pdf when x = is sent is given by /o(/) = J2jLiPj^ {j ~ ^"-"S^j ^iid by 
/i(0 = Sjii Qj^ ~ log j'^ when x = 1 is sent over the channel. 

Now, assuming that when yj is observed at the channel output, the approximate LLR is 
calculated by Ij = g{yj) = aj, the conditional pdfs of / are: 



/o(0 



— a. 



M 



Inserting ([HI) and ([T5]) in (H)) gives 



1 



(14) 
(15) 

(16) 



Taking |^ reveals that aj = log ^ maximizes C*j for all 1 < j < Af since < for all 
1 < i < M and = for all 1 < j < M and 1 < k < M and j 7^ k. These values of 

a/s are equal to the true LLR of (fT3l) . Thus, the maximizing point is only given by true LLRs. 
Noticing that these results are valid for each equivalent bit-channel i of the BICM and since 
max C = Xlili Hiax^. Ci in ([5]), the theorem is proved. 

Appendix B 
Proof of Theorem [2] 



Denote = E^^;:^.[U'^\^] and y^*^ = E^^;^^[y\x] for b e {0, 1}. Then A can be written 



as 



2 ^0 
--1 - -E (0 

2 yo 



log2(l + 
logs ( 1 + e 



log2(l + e^i 

log2 ( 1 + e 



2 yi 



By using ([7]) and with some abuse of notation we write 

1 iv(*) 



b=0 k=l 



log2 (l + e(-i)'"^«4^''y^)+^"')) 



For each fixed y5 G C^*\ it is clear that is a linear function of a[.*^ and Noticing that 
the function logg (1 + exp(-)) is convex and twice differentiable, it can be deduced that log2(l + 
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exp((— l)''+^((aj!\ Yb) + 6^'^))) is also a convex function of a^!^ and h^^ . The convexity is also 
preserved under expectation. Thus, E. (i)Jlog2(l + exp ((— l)''+^((a^'-', yt) + 61*^)))] is also 
convex which makes Q concave with respect to a[!^ and l)^' for all A; = 1, ... , A^*^*). 

References 

[1] S. Lin and D. J. Costello, Error Control Coding, Second Edition. Upper Saddle River, NJ, USA: Prentice-Hall, Inc., 2004. 
[2] C. Berrou, A. Glavieux, and R Thitimajshima, "Near Shannon limit error-correcting coding and decoding: Turbo-codes (1)," 

in Conf. Rec. IEEE ICC '93, vol. 2, Geneva, Switzerland, May 1993, pp. 1064-1070. 
[3] R. G. Gallager, "Low-Density Parity-Check codes," Ph.D. dissertation, MIT press, Cambridge, MA, 1963. 
[4] I. Land, P. A. Hoeher, and U. Sorger, "Log-likelihood values and Monte Carlo simulation: some fundamental results," in 

Proc. Int. Symp. on Turbo Codes and Related Topics, Brest, France, Sept 2000. 
[5] T. J. Richardson, M. A. Shokrollahi, and R. L. Urbanke, "Design of capacity-approaching irregular low-density parity-check 

codes," IEEE Trans. Inf. Theory, vol. 47, no. 2, pp. 619-637, Feb. 2001. 
[6] A. Abedi and A. K. Khandani, "Invariance properties of binary linear codes over a memoryless channel with discrete input," 

IEEE Trans. Inf. Theory, vol. 53, no. 3, pp. 1215-I2I8, Mar 2007. 
[7] J. Hagenauer, "Viterbi decoding of convolutional codes for fading- and burst-channels," in Zurich seminar on digital 

communications, 1980. 

[8] J. Hou, P. H. Siegel, and L. B. Milstein, "Performance analysis and code optimization of low-density parity-check codes on 

Rayleigh fading channels," IEEE J. Sel. Areas Commun., vol. 19, no. 5, pp. 924-934, May 2001. 
[9] F. Tosato and P. Bisaglia, "SimpUfled soft-output demapper for binary interleaved COFDM with application to HIPERLAN/2," 

in Proc. of the 2002 IEEE International Conference on Communications (ICC'02), 2002, pp. 664—668. 
[10] J. K. Kwon, S. Park, and D. K. Sung, "Log-likelihood ratio (LLR) conversion schemes in orthogonal code hopping 

multiplexing," IEEE Commun. Lett., vol. 7, no. 3, pp. 104-106, Mar. 2003. 
[11] ETSI TS 101 475, "Broadband Radio Access Networks (BRAN); HIPERLAN Type 2; Physical (PHY) layer, vl.2.2," 2001. 
[12] R. Yazdani and M. Ardakani, "Linear LLR approximation for iterative decoding on wireless channels," IEEE Trans. Commun., 

vol. 57, no. 11, pp. 3278-3287, Nov. 2009. 
[13] G. Caire, G. Taricco, and E. Biglieri, "Bit-interleaved coded modulation." IEEE Trans. Inf. Theory, vol. 44, no. 3, pp. 

927-946, 1998. 

[14] S.-Y. Chung, G. D. Forney, T. J. Richardson, and R. Urbanke, "On the design of low-density perity-check codes within 
0.0045 dB of the Shannon limit," IEEE Commun. Lett., vol. 5, no. 2, pp. 58-60, Feb. 2001. 

[15] M. Luby, M. Mitzenmacher, A. ShokroUahi, and D. Spielman, "Efficient erasure correcting codes," IEEE Trans. Inf. Theory, 
vol. 47, pp. 569-584, 2001. 

[16] J. Hou, P. H. Siegel, L. B. Milstein, and H. D. Pflster, "Capacity-approaching bandwidth-efficient coded modulation schemes 
based on low-density parity-check codes," IEEE Trans. Inf. Theory, vol. 49, no. 9, pp. 2141-2155, Sept. 2003. 

[17] A. Sanaei and M. Ardakani, "LDPC code design considerations for non-uniform channels," IEEE Trans. Commun., Jan. 
2010. 

[18] M. G. Luby, M. Mitzenmacher, M. A. ShokroUahi, and D. A. Spielman, "Improved low-density parity-check codes using 
irregular graphs," IEEE Trans. Inf Theory, vol. 47, no. 2, pp. 585-598, Feb. 2001. 



11 



8-AM 


SNR 


Bit 1 


Bit 2 


Bit 3 


7.88 dB 


a<^' = 1.328 


^ = 0.612 


= 2.046 


^ = 0.328 


bf^ = 2.273 


a^^' = -0.482 


4=^) = -0.909 


cp^ = -3.928 


21.03 dB 


a^^' = 8.538 


^ = 0.825 


' = 3.098 


' = 0.384 


b'f = 2.513 


= -1.528 


6(3) = _2.357 


= -2.547 


16-QAM 


SNR 


Bit 1 


Bit 2 


5.02 dB 


= 1.262 


' = (0.868, -0.200) 


6^ ' = -1.257 



TABLE I 

Optimized piecewise linear LLR parameters at different SNRs for 8- AM and 16-QAM when K = 0. 
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Fig. 1. The 8-AM and 16-QAM constellation points and Gray mapping. 
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(c) i = 3 

Fig. 2. True bit LLR values LS^^ as functions of the channel output y for the 8-AM at SNR= 7.88 dB. Also, the piecewise 
linear LLR approximations with optimized parameters of Table U are depicted for K — 0. 
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Fig. 3. Comparison between the BER of a randomly constructed (3, 4)-regular LDPC code of length 15000 decoded by true 
and approximate LLRs on the Rayleigh fading channel (K = 0). The approximate LLR parameters are reported in Table U 



